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Abstract. 

We report on the analysis of a BeppoSAX observation of IRAS 13197-1627. The complexity of the broad band 
(0.3-100 keV) spectrum of the source prevents us from selecting from two possible best-fit models, i.e. one based 
on a partial covering and the other on a reflection dominated scenario. Whatever the best fit model, the data 
imply that: 1) the primary spectrum is heavily obscured by at least two absorbers with column densities of the 
order of 5x10'^'' cm~^ and >2xl0^* cm~'^; 2) the intrinsic 2-10 keV luminosity of the source is at least ~2xl0'*'' 
erg s~^ thus making IRAS 13197-1627 a type-II QSO (more precisely a type 1.8 QSO) rather than a Seyfert 
1.8 galaxy as previously classified. Furthermore, there is marginal evidence of an absorption line at E ~7.5 keV 
suggesting a possible outflow with t;~0.1c similar to recent findings on other QSOs and Seyferts. 
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1. Introduction 

Much of the research on Active Galactic Nuclei (AGNs) 
in recent years emphasized the search for Type-II QSOs. 
Unified models of AGNs (e.g. Antonucci 1993) explain 
the observed differences between broad (Type-I) and nar- 
row (Type-II) emission-line AGNs as obscuration and 
viewing angle effects rather than intrinsic physical dif- 
ferences. Although rather successful in describing nearby 
Seyfert galaxies, unified models seem to fail when ex- 
tended to their higher-luminosity "relatives" , the QSOs 
(hereinafter, we shall broadly define a Seyfert as having 
L(2-iokcV) <10^'^ergs-^ andaQSOifL(2_iokoV) >10''''erg 

On the basis of unified models, the existence of a large 
population of type-II QSOs, or broadly speaking of a large 
population of heavily-absorbed high-luminosity AGNs has 
to be expected. However, because of the intrinsic difficul- 
ties to find them, only a limited number of such objects has 
been identified. Their rarity is somehow embarrassing con- 
sidering that they are supposed to contribute the most to 
the total emission of the cosmic X-ray background. This is 
now supported both observationally from recent Chandra 
and XMM-Newton deep-fields and theoretically from syn- 
thesis models of the cosmic X-ray background (see e.g. the 
recent review by Fabian 2003). 

In this paper we report results from the BeppoSAX 
narrow field instruments (Boella et al. 1997a) observa- 
tion of IRAS 13197-1627. According to these, we find that 
the source can be identified as the brightest Type-II QSO 
known to date. Given the above argument, it is thus an in- 



teresting case-study to understand the detailed absorption 
properties of this type of objects. 



Table 1. Basic parameters of IRAS 13197-1627. (1) de 
Robertis et al. 1988; (2) Aguero et al. 1994; (3) Theureau 
et al. 1998. 



Hubble Type 


S? (1) 


Seyfert Type 


1.8 (2) 


Red shift 


0.01654 (3) 


r.a. (2000) 


13 22 24.462 


Dec. (2000) 


-16 43 42.91 




14.46 



IRAS 13197-1627, also named MCG-3-34-64i, has been 
firstly classified as a Seyfert 2 galaxy ( Osterbrock & de 
Robertis 1985, de Robertis et al. 1988) and, afterward, as 
a type 1.8 Seyfert by Aguero et al. (1994) and Young et 
al. (1996). Its basic parameters are given in Table 1. 



The source has been often misidentified with MCG-3-34- 
63, an edge-on galaxy -2' north-est of IRAS 13197-1627. We 
checked that the X-ray emission is indeed coming from MCG-3- 
34-64 while MCG-3-34-63, is not detected in 0.1-10 keV band. 
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IRAS 13197-1627 

Bin time; 4200. 




2x10* 4x10* 6x10* 8x10* 10^ 

Time (s) 



Fig. 1. MECS 2-5 keV (upper panel) 5-10 keV (middle panel) 
light curves and their ratio as a function of time (lower panel) . 

The morphological classification of its host galaxy is 
unclear. In optical images it is extended over ~30"x30" 
with only a few indications, if any. of spiral structures (dc 
Robertis et al. 1988). A broad component in the Ha line 
has been detected in polarized light by Young et al. (1996). 

IRAS 13197-1627 was first observed in X-rays by 
ASCA (Ueno 1995). The ASCA spectrum was extremely 
steep (the photon index was r~3.0±0.3), with a high ab- 
sorption column density of Nh'^4x10^^ cm~^ and with 
a soft fainter component emerging at energies below ~3 
keV. ASCA also detected an FeKa line at E~6.4 keV with 
an equivalent width EW--350 eV (Ueno 1995). 

The present BeppoSAX dataset has been previously 
analyzed by Risaliti (2002) in a statistical study of a sam- 
ple of 20 bright type-II Seyfert galaxies. However, Risaliti 
(2002) used only a predefined set of 3 spectral models and 
did not attempt a detailed modeling of all the spectral 
complexities found in this particular source. The reduced 

of the best-fit reported by Risaliti (2002) was ~1.3, 
much larger than those reported here. The present work 
thus is complementary to and expands Risaliti's results. 

2. Observation and data reduction 

BeppoSAX observed IRAS 13197-1627 on July 22, 1998. 
The LEGS (Parmar et al. 1997), MECS (Boella et al. 
1997b) and PDS (Frontera et al. 1997) exposure times 
were of 14 ks, 44 ks and 19 ks respectively. 

Data reduction has been performed using the standard 
mission specific software SAXDAS based on the FTOOLS 
package. 

Two sources are visible at about 20' from the target 
in the MECS 2-10 keV image. None is detected by the 
LECS so we can exclude that they contaminate the IRAS 
13197-1627 counts at low energies. These sources are also 
not detected by the MECS in the 5-10 keV energy range. 
For this reason they should not contribute significantly to 
the PDS counts. We checked in public archived databases 




1 10 
Energy (keV) 



Fig. 2. Data-to-modcl ratios obtained when fitting the ~0.3- 
100 keV data with a power-law model with Galactic absorp- 
tion. 

that no other sources with known strong X-ray (2-10 keV) 
emission are present in the sky field covered by the PDS 
instrument (1.3° FWHM, Frontera et al. 1997). 

Source counts have been extracted from circular re- 
gions centered on IRAS 13197-1627. The LECS and MECS 
extraction regions have radii of 8' and 4' respectively. The 
LECS, MECS and PDS spectra have been rebinned so as 
to sample the energy resolution of the instruments using 
grouping files produced by the ASI Science Data Center 
(ASDC) and taking care of having at least 20 counts for 
each obtained energy bin. To subtract the background, we 
used standard PHA files accumulated from observations 
of empty sky regions and produced by the BeppoSAX 
team. We also checked that the results obtained using lo- 
cal backgrounds are consistent with the ones presented 
here. Most recent calibration files have been used for the 
spectral analysis. 

Fig. 1 shows the MECS 2-5 keV (upper panel) and 
5-10 keV (middle panel) lightcurves. The source flux is 
found to be roughly constant during the first half of the 
observation, while it is significantly (at a confidence level 
~95% ) variable by up to 40% on timescales down to 5-10 
ks during the second half. These variations are observed 
both in the 2-5 keV band where the source spectrum is 
strongly affected by absorption and in the 5-10 keV energy 
range where the source should be seen almost directly. 

As implied by the hardness ratios shown in Fig. 1 
(lower panel), we found no significant spectral variations 
(with a confidence level >99%) during the observation. 
We thus performed a time averaged spectral analysis col- 
lecting events from the whole period. 

3. Spectral analysis 

As shown by the data-to-model ratios plotted in Fig. 2, 
the broad-band (~0.3-100 keV) spectrum of IRAS 13197- 
1627, once fitted with a simple power-law, clearly shows 
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Table 2. SME model. (1) model name; (2) photon index; (3) absorption column density; (4) energy of the absorption 
edge; (5) optical depth of the absorption edge; (6) x^/degrees of freedom. 



Model 



(1) 



SME 



(2) 



Nh 
10^^ cm" 
(3) 



keV 
(4) 



Tedge X^/d.O.f. 



(5) 



^+0.17 
"'-0.15 



50.9+^;^ 



7 00+°-^^ 



(6) 



Energy (keV) 



1 ' 1 ' 1 ' 1 ' 1 
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Fig. 3. The 3-10 koV residuals (expressed in terms of standard Fig. 4. Confidence contour levels (68, 90 and 99%) for the edge 
deviation) for the SM. For clarity, data have been rebinned so optical depth t versus its energy centroid. 
as to have a 3a significance for each plotted data point. 



three main spectral features: 1) a soft excess below ~2- 

3 kcV; 2) a low-cncrgy cut-ofF due to strong absorption 
between ~3-5 keV; 3) a sharp and deep drop of the counts 
at E~7 keV. 

3.1. Simple Spectral Parametrization 

To account for the strong absorption and the soft excess 
we added to the simple power-law model a cold absorber 
(Nh~5±1x10^^ cm~^) plus a power-law emerging at low 
energies. In a classical Seyfert 2 spectrum the latter com- 
ponent would be interpreted as due to scattering of the pri- 
mary emission (hereinafter we will mention to this model 
as the "scattering model" SM). The data-to-modcl resid- 
uals (Fig. 3) for SM still show remaining complexities in 
the ^4-8 keV energy range: an emission feature at 6.4 
keV, a deep absorption edge at ^7 keV and possible other 
components arising between ~4-6 keV. 

To better model the data, we first tried to obtain a sim- 
ple parametrization of all the features disregarding their 
physical implication. 

The quality of the fit increases drastically after the 
addition of an absorption edge to the SM (Ax^=24.3 
for 2 additional parameters, see SME model in Table 2). 
According to the F-test, this absorption edge is significant 




^ I I I I I I 

5 10 
Energy (keV) 

Fig. 5. Residuals in the 3-10 keV energy range for the SME 
model (data binned as in Fig. 3) . 

at >99.99% confidence level (see Fig. 4). The energy of the 
edge (E~7.0±0.2 keV) rules out (at >99%) the presence 
of iron ionized at levels larger than FeX. If produced by 
the measured absorption column (Nh~5x10^'^ cm^^) the 
depth of the edge would require an iron overabundance 
ranging between 4 and 8 times the solar value. 
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Table 3. Partial Covering Model. (1) model name; (2) photon index; (3) column density of the whole absorber; (4) 
covering fraction of the partial absorber; (5) column density of the partial absorber; (6) energy of the FeKa emission 
line; (7) equivalent width (EW) of the FeKa line; (8) x^/degrees of freedom. 



Model 
(1) 


r 

(2) 


Nn.whole 

10^^ cm-^ 
(3) 


Cf 
(4) 


Nh, partial 

(5) 


EpeKa 

keV 
(6) 


EW 
eV 
(7) 


PC 


, 7+0.34 


47.8i|;i 


94+° °i 

"■^^ — 0.03 




6 42+°'^^ 

"•^•'-0.12 


228i;?« 



xVd.o.f. 
(8) 



82.9/90 



Moreover, the data-to-modcl residuals for the SME 
model still show significant features between ~4 and 8 keV 
(see Fig. 5). These features could be fitted by Gaussian 
emission lines (one broad at £^^4.5 kcV and one narrow 
at E~6.4 keV) and/or by Gaussian absorption fines (at 
E~6 keV and E~7.5 keV), afi of which significant at more 
than 90% confidence. 

In summary, the "simple spectral parametrization" 
demonstrates that a single absorber cannot describe weU 
the data above -^2 keV. 

3.2. Partial Covering 

A possible explanation for the deep edge detected at ~7 
keV could be the occurrence of an additional heavy ab- 
sorption that covers only partially the primary emission. 
This model fits well the data (model PC in Table 3 and 
Fig. 6) and the parameter values obtained indicate that 
about 95% of the primary emission is absorbed by cold 
matter with Nh~2x10^^ cm~^. The photon index of the 
continuum is quite steep (r~3.2±0.3) while the FeKc, line 
is detected (at a confidence level >99%) at E^6.4±0.1 
keV with EW~228±180 eV. With this model the lu- 
minosity of the source is found to be L2_iokeV~2xlO^^ 
erg s^^. The observed FeK^ line EW is lower than ex- 
pected for the observed column densities (an absorbing 
column of Nh~2x10^* cm~^ with a ~95% coverage plus 
a Nh^SxIO^'' cm~^ with 100% coverage should imply an 
EW~700 eV, Makishima et al. 1986). 

3.3. Pure Reflection 

Alternatively to the PC model, the edge at 7 keV could 
be produced by a reflection component, which dominates 
the overall spectral shape. 

If we add to any of the above mentioned models a cold 
reflection component (namely PEXRAV model in XSPEC, 
Magdziarz & Zdziarski 1995), we obtain that the R param- 
eter that define the relative importance of the reflected 
continuum with respect to the direct one, reaches very 
high values (i?>10), thus reproducing a completely re- 
flected spectrum. 

We thus tried a model where the primary continuum 
is completely hidden for us and we only see the contin- 
uum reflected from the disk plus the associated FeKa line. 



Energy (keV) 

Fig. 6. Data-to-model residuals for the PC model (data binned 
as in Fig. 3). 



Energy (keV) 

Fig. 7. Data-to-model residuals for the PR model (data binned 
as in Fig. 3). 



This latter component has been parametrized using a rel- 
ativistic model (the Laor model in XSPEC, see Fabian 
et al. 2002 for a recent review on this topic). Since the 
data clearly show a low energy cut-off, we also added an 
absorption column to the model plus the associated nar- 
row Gaussian emission line. To prevent the parameters 
of this model to diverge we fixed some values of the re- 
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Table 4. Pure Reflection model. (1) model name; (2) photon index; (3) absorption column density; (4) cut-off energy; 
(5) energy of the relativistic FeKa line; (6) external radius of the region of emission of the relativistic line expressed in 
units of gravitational radii; (7) equivalent width of the relativistic line; (8) energy of the narrow FeKa line; (9) equivalent 
width of the narrow line; (10) /degrees of freedom. ^ pegged parameter. Fixed parameter. 

Model r Nh Ecut-off ELaor Tout EWLaor E„arr. EW„arr. X^/d.O.f. 

10^^ cm-^ keV keV keV keV keV eV 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



PR 2.30t°;|» 41.2iiJJ rettf'" 6.40f 6.05l«jJ S.Slifg 6.40l°;i° 356^111 66.8/88 



flection component: metal abundance, R and inclination 
were fixed to solar, 1 and 60° respectively. Tlie only free 
parameters of the continuum were thus the power-law pho- 
ton index, the high energy cut-off and the normalization. 
Similarly, we fixed many parameters of the relativistic line 
as follow: line energy E=6.4 keV, emissivity index i=-3, 
inclination anglc=60°, internal radius of the emitting re- 
gion=1.25 gravitational radii r^. The only free parameters 
of the Laor model thus were the external radius of the line 
emitting region (rout) and normalization. 

This model well describes the data (PR model in Table 
4 and Fig. 7) and the inferred intrinsic luminosity of the 
source is found to be L2_iokoV>2x lO**^ erg s~^. The fit 
results do not depend on the energy of the relativistic line. 
In fact, values of ELaor from 6.4 to 6.9 keV give fits of the 
same quality (Ax^< 0.2). The EW of the relativistic hne 
is rather large (EW~3.8 keV) and it is in agreement with 
what is expected from a pure reflection spectrum (Matt 
et al. 1996). The outer radius of the line emitting region is 
found to be ~5-7 r^, a rather extreme value. The energy 
and EW of the narrow line are in perfect agreement with 
what expected for a lino produced via transmission in a 
cold absorber with the measured Nh of ~4xl0^^ cm~^ 
(Makishima 1986). 

3.4. An absorption line at 7.5 keV? 

In the data-to-model residuals both for the PC and PR 

model there is a marginal evidence of an additional ab- 
sorption feature at E~7.5 keV. We checked this hypothesis 
by adding an absorption Gaussian line to the two models. 

The quality of the fit drastically increases in the case of 
the partial covering (Ax^=15.2 for 2 additional parame- 
ters, Fig. 8) while the feature seems to be marginal adopt- 
ing the PR model as baseline (Ax^=5.4). 

The feature is found to be at E~7.5±0.2 keV and to be 
narrow {a< 0.47 keV) while its EW is ~-330±180 eV Its 
narrowness, in particular, seems to exclude that it could be 
associated to an unresolved additional edge due to warm 
matter. 

This absorption line could be associated with resonant 
absorption due to He-like and/or H-like Fe. We tested this 
hypothesis modeling the feature with the SIABS model 




k7 7.2 lA 7.6 7.B 



Energy (keV) 

Fig. 8. Confidence contour levels for the absorption line at 
~7.5 keV in the PC model. 

(Kinkhabwala et al. 2003) assuming it is due to He- like Fc. 
We obtained a good quality fit with physically plausible 
parameters, i.e. NFe~10^^ cm~^, outward velocity t;'~0.11c 
and dispersion At;'~4500 km s~^. 

4. Discussion 

The overall 0.3-100 keV spectrum of IRAS 1319-1627 is 
dominated by cold absorption and by an extra soft com- 
ponent that rises below ~3 keV. This soft component is 
well modeled by a power-law with the same photon index 
of the hard component. Thus, as a first approximation, 
this spectrum appears as a very common (i.e. highly ab- 
sorbed) X-ray spectrum of a Seyfert 2 galaxy where the 
low energy band is dominated by a scattered component 
(Smith and Done 1996, Turner ct al. 1997, Risahti 2002). 

However, unlike other spectra of known Seyfert 2s, 
there is a major unexpected feature in the spectrum of 
IRAS 13197-1627 that is a very deep absorption edge at 
-^7 keV. The edge energy is consistent with it being pro- 
duced in a cold medium and rules out the presence of any 
warm material. The depth of the edge requires an absorber 
with Nh~2x10^^ cm~^, i.e. ^3-4 times larger than what 
found by modeling the low energy cut-off with a single ab- 
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sorption column. Thus, to explain this edge, we are forced 
to invoke either that the source is partially covered or that 
it is completely obscured and its spectrum is dominated 
by the reflected component. 

4.1. The Pure Reflection and Partial Covering 
Scenarios 

The pure reflection scenario seems rather extreme since 
a strong relativistic line (EW^^ 3.8 keV) produced within 
~7 gravitational radii is required to fit the data. Such an 
extreme line recalls the one that has been recently invoked 
for the narrow line Seyfert 1 MCG-6-30-15 (Wilms et al. 
2001). 

Bassani et al. (1999) suggested to use the ratio Rct 
between the absorption corrected OIII[5007A] flux and the 
observed 2-10 kcV flux to identify Compton-thick sources. 
In their scheme, Compton-thick sources have Rct<0.1. 
In the present case the 2-10 keV flux is ~4xl0~^^ erg 
cm^^ and the absorption corrected F[5077A] flux 
is ~4xl0~^^ erg cm~^ s~^ (Young et al. 1996) yield- 
ing RcT~l. Also using the ASCA flux (Ueno 1995) one 
obtains Rct"^0.5. Therefore, according to this critcriuni, 
IRAS 13197-1627 should not be a Compton-thick source, 
making the pure reflection scenario even more implausible. 

The partial covering scenario appears to be more likely 
not last for its simplicity. Nevertheless also in this case 
there are difficulties in interpreting the best fit parameters. 
First of all, the best fit photon index r'~3 is rather steep 
for Seyfert galaxies (Smith and Done 1996, Turner et al. 
1997) and it is only marginally consistent with what found 
for the most extreme Narrow Lines Seyfert 1 (NLSYl, 
Boiler Brandt & Fink 1996, Vaughan et al. 1999), i.e. the 
radio quiet ACNs with the known steeper spectra. It is 
worth noting here that another possible similitude with 
this class of sources is the observed variability (^40%) 
down to timescales of ~5-10 ks. Another inconsistency in 
the partial covering scenario is that the EW of the FeKa 
line is lower than what expected for the observed absorp- 
tion (see Sect. 3.2). 

Present results show interesting analogies with what 
found for NLSYls. This class of objects are generally char- 
acterized in X-rays by a strong soft excess, a steep power- 
law and extreme variability (Boiler et al. 1996, Brandt 
et al. 1997, Brandt & Boiler 1998, Vaughan et al. 1999). 
Recent XMM-Newton results have added to this list, a 
new "typical" feature: a sharp spectral drop above 7 keV 
(Boiler et al 2001, Pounds et al 2003a, Boiler et al. 2003). 
Apart from the soft X-rays (that are here absorbed), IRAS 
13197-1627 shows all the above mentioned characteristics 
of NLSYls. The spectral complexity measured in IRAS 
13197-1627 is indeed very similar to what recently ob- 
served in MKN 766 (Pounds et al. 2003a), IH 0707-495 
(Boiler et al. 2002) and IRAS 13224-38098 (Boiler et al. 
2003a) where also in these cases the authors were left with 
two possible scenarios that described the data equally well: 



either a refiection dominated or a partial covering sce- 
nario. 

A theoretical explanation for the reflection dominated 

scenario has been recently proposed by Fabian et al. 
(2002). The main assumption of this model is that, in 
high accretion rate systems like NLSYls, instabilities in 
the accretion may lead to the formation of multiple reflec- 
tors in the form of cold, dense material clumped into deep 
sheets or rings. 

Within the partial covering scenario, it could be that, 
as proposed by Pounds et al. (2003a) for MKN 766, 
dense ionized clouds or ejecta close to the super-massive 
black hole are responsible for the deep edge observed with 
XMM-Newton. Also in this case our data are in good 
agreement with the model with the exception of the ion- 
ization state of the partially absorbing matter, our data 
requiring near-neutral Fe. 

To disentangle between these two hypothesis better 
quality data above 3 keV (e.g. from XMM-Newton) are 
needed. Nevertheless, it is worth noting here that both 
scenarios imply that the intrinsic 2-10 keV luminosity of 
the source is > 2x10** erg s'^ This makes IRAS 13197- 
1627 the nearest and brightest type-1.8 QSO known to- 
date. This source thus offers a unique possibility to study 
the spectral properties of this class of objects and to test 
unified models at the bright end of the luminosity function 
of Seyfert galaxies. This is of particular importance given 
that these objects are thought to contribute the most of 
the cosmic X-ray background. IRAS 13197-1627 is thus an 
important case study to understand what the absorption 
characteristics of luminous QSO could be. 

4.2. Possible Origin of the Complex Absorption 

As stated above, the spectral analysis shows that the 
source has heavy and complex absorption that cannot be 
modeled by a single absorber. 

We speculate that the most plausible scenario is the 
one in which the lower column is due to the putative dusty 
torus while the heavier absorption is due to dense clouds, 
blobs or winds close to the X-ray emitting region. The 
presence of a dusty torus, in fact, would explain the differ- 
ent columns measured in X-rays (Nh'~5x10^^ cm~^) and 
in optical (the Eb-v in the direction of the nucleus implies 
aNH~1.5xlO^^ cm"2, Lumsden et al. 2001). On the other 
hand, dense blobs or clouds close to the X-ray source may 
either partially cover or block the nuclear emission with- 
out obscuring the component reflected by the accretion 
disk. 

It is worth noting here that this picture would also 
be able to explain, if confirmed, the absorption feature at 
~7.5 keV. This, in fact, could be interpreted as He- or 
H-like Fe resonant absorption by outflowing matter. The 
inferred velocity of this outflow then should be ~0.1c, sim- 
ilar to what found in PG1211-M43 (Pounds et al. 2003). 
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5. Conclusion 

The BeppoSAX X-ray spectrum of IRAS 13197-1627 is 
so complex that we could not find a single best-fit model 
and were left with two possible scenarios: i) one where the 
primary emission is completely blocked and the resulting 
spectrum is completely refiection dominated and ii) one 
where the spectrum is due to a partial covering of the 
X-ray source. 

Despite this ambiguity, a firm conclusion of this work 
is that any modeling of the X-ray spectrum implies an 
intrinsic luminosity of at least L2_iokeV~2.0x 10^^ erg s~^ 
thus making IRAS 13197-1627 the nearest and brightest 
example of a type-II QSO. 

Moreover, both scenarios require a complex absorption 
structure composed of at least two heavy absorbers. We 
speculate that these could be a dusty torus and clouds, 
blobs or winds closer to the nucleus. Finally, if the ab- 
sorption line marginally detected at 7.5 keV is identified 
with resonant absorption from He- or H-likc Fc, then we 
infer that some of matter (e.g. the clouds or the wind) 
should have an outward velocity of the order of ~0.1c. 
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